Anomalous Excitation Spectra of Frustrated Quantum Antiferromagnets 
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We use series expansions to study the excitation spectra of spin- 1/2 antiferromagnets on 
anisotropic triangular lattices. For the isotropic triangular lattice model (TLM) the high-energy 
spectra show several anomalous features that differ strongly from linear spin- wave theory (LSWT). 
Even in the Neel phase, the deviations from LSWT increase sharply with frustration, leading to 
roton-like minima at special wavevectors. We argue that these results can be interpreted naturally 
in a spinon language, and provide an explanation for the previously observed anomalous finite- 
temperature properties of the TLM. In the coupled-chains limit, quantum renormalizations strongly 
enhance the one-dimensionality of the spectra, in agreement with experiments on CS2CUCI4. 

PACS numbers: 75.10.Jm 



One of the central problems in quantum magnetism 
is understanding the properties of two-dimensional (2D) 
spin-1/2 Heisenberg antiferromagnets (HAFM's). A 
question of particular interest is whether the interplay 
between quantum fluctuations and geometrical frustra- 
tion can lead to unconventional ground states and/or 
excitations. Candidate materials which have recently 
attracted much attention include CS2CUCI4 Q and k- 
(BEDT-TTF) 2 Cu 2 (CN) 3 0. 

If the ground state is magnetically ordered, the sys- 
tem must have gapless magnon excitations, which at suf- 
ficiently low energies are expected to be well described 
by semiclassical (i.e., large- S) approaches like spin- wave 
theory (SWT) and the nonlinear sigma model (NLSM). 
However, if the magnon dispersion at higher energies de- 
viates significantly from the semiclassical predictions, it 
is possible that the proper description of the excitations, 
valid at all energies, is in terms of pairs of S = 1/2 
spinons. In this unconventional scenario the magnon is 
a bound state of two spinons, lying below the two-spinon 
(particle-hole) continuum. 

In this Letter we use series expansions to calculate the 
magnon dispersion of 2D frustrated 5 = 1/2 HAFM's. 
Our main finding is that for the triangular lattice model 
(TLM) the dispersion shows major deviations from lin- 
ear SWT (LSWT) at high energies (Fig. P). We argue 
that these deviations can be qualitatively understood in 
terms of a two-spinon picture Q, provided the spinon 
dispersion has minima at Ki/2, where K$ is a magnetic 
Bragg vector. Based on this interpretation of the TLM 
spectra, we propose an explanation for the anomalous 
finite-temperature behavior found in high temperature 
series expansion studies 0. 

Both qualitatively and quantitatively the deviations 
from LSWT found here for the TLM are much more pro- 
nounced than those previously reported [E El 0,0] f° r trie 
high-energy spectra of the square lattice model (SLM). 



We point out that the deviations from SWT increase in 
the Neel phase too, upon adding frustration to the SLM. 
We further consider the limit of our model relevant to 
CS2CUCI4, and show that the calculated excitation spec- 
tra are in good agreement with experiments 0. 

Model— We consider a S = 1/2 HAFM on an 
anisotropic triangular lattice, with exchange couplings 
Ji and J 2 (Fig. Efa)). This model interpolates between 
the SLM (Ji = 0), TLM (Ji = J 2 ), and decoupled 
chains (J 2 = 0). Classically the model has Neel order for 
•h _ ^2/2 with q — 7T, and helical order for Ji > J 2 /2 
with q — arccos (— J 2 /2 Ji), where q (2q) is the angle be- 
tween nearest-neigbor spins along J 2 (</])• The phase 
diagram for S = 1/2 was studied in Rcfs. 103 

Series expansion method. — In order to develop series 
expansions for the helical phase (see Ref . for expansions 
for the Neel phase spectra) , we assume that the spins or- 
der in the xz plane, with nearest-neighbor angles q and 
2q as defined above; q (generally different from the classi- 
cal result) is determined by minimizing the ground state 
energy. We now rotate all the spins so as to make the 
ordered state ferromagnetic, and introduce an anisotropy 
parameter in the Hamiltonian H(\) = Hq + XV, so that 
H(0) is a ferromagnetic Ising model and H (1) is the spin- 
rotation invariant Heisenberg model @. We use linked- 
clustcr methods to develop series expansions in powers of 
A for ground state properties and the triplet excitation 
spectra. The calculation of the spectra is particularly 
challenging as S z is not conserved. Due to single-spin 
flip terms in V, the one-magnon state and the ground 
state belong to the same sector, and the linked-cluster 
expansion with the traditional similarity transformation 
|11| fails. To get a successful linked-cluster expansion, 
one has to use a multi-block orthogonality transforma- 
tion We have computed the series for ground state 
properties to order A 11 , and for the spectra to order A 9 
for J 1 = J 2 and to order A 8 otherwise. The properties 
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for A = 1, discussed in the following, are obtained from 
standard series extrapolation methods. 
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FIG. 1: Excitation spectrum for the TLM (Ji = J2) along the 
path ABCOAQD shown in Fig. 0] The high-energy spectrum 
is strongly renormalized downwards compared to the LSWT 
prediction (full line). Note the "roton" minima at B and D 
and the flat dispersion in the middle parts of CO and OQ. 

Square lattice model. — Several numerical studies 0, 
IE Hi have reported deviations from SWT for high- 
energy excitations in the SLM. While 1st order SWT 
(i.e., LSWT) and 2nd order SWT predict no dispersion 
along (tt — x, x), 3rd order SWT finds a weak dispersion, 
with the energy at (tt, 0) ~ 2% lower than at (tt/2,tt/2) 
d,E3,0|. In contrast, the most recent Quantum Monte 
Carlo (QMC) and series expansion [8( studies find 
the energy at (tt, 0) to be ~ 9% lower. These deviations 
from SWT have been interpreted 0, 0, in terms of 
a resonating valence-bond (RVB) picture, in which the 
ground state is described as a 7r-ffux phase modi- 
fied by correlations producing long-range Neel order , 
and the Goldstone modes (magnons) are bound states of 
a particle and a hole spinon |l5|, [l6j . The magnon dis- 
persion, calculated using a random phase approximation 
(RPA), has local minima at (tt, 0) 0,0] (in qualitative 
agreement with the series/ QMC results), the locations of 
which are intimately related to the fact that the spinon 
dispersion has minima at (7r/2,7r/2) [l5llr^ |. 

Frustrated square lattice model. — The Neel phase per- 
sists up to J\j Ji = {J\/J2)c ^ 0.7, after which the 
system enters a dimerized phase [$]. As the frustra- 
tion J\j Ji is increased towards (J\j J-i) Cl the local mini- 
mum at (tt, 0) becomes more pronounced (see Fig. |3J); for 
J1/J2 — 0.7 the energy difference between (tt/2, it/2) and 
(tt, 0) has increased to ~ 31%. In contrast, LSWT pre- 
dicts no energy difference 0] . These results lend further 
support to the RVB/flux-phase picture. The locations of 
the Bragg vectors, roton minima, and spinon minima in 
the Neel phase are shown in Fig. [2Jb). 
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FIG. 2: (a) Exchange constants J\ and J2 for the 5 = 1/2 
HAFM on the anisotropic triangular lattice. The model can 
also be viewed as a square lattice with an extra exchange along 
one diagonal, (b) Brillouin zone for the frustrated SLM, in 
standard square lattice notation (the frustrating Ji bonds are 
taken to lie along the +45° directions in real space) . The exci- 
tation spectra in Fig. [3] are plotted along the bold path. Also 
shown are the locations of the Bragg vectors (filled squares) 
and the local "roton" minima (filled circles) in the 5 = 1 
dispersion, and the global minima of the 5=1/2 spinon dis- 
persion (open squares) in the Neel phase. Note that the 90° 
rotation invariance present for Ji = is lost for J\ > 0. 
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FIG. 3: Excitation spectra in the Neel phase. As J\jJ% 
is increased, the local "roton" minimum at (tt, 0) becomes 
more pronounced, and the energy difference between (tt, 0) 
and (7r/2,7r/2), which LSWT predicts to be zero, increases. 



Triangular lattice model. — Next, we consider the TLM 
(J\ = J2), whose ground state has 120° ordering be- 
tween neighboring spins 0] . Fig. ^ shows the excita- 
tion spectrum plotted along the path ABCOQD in Fig. 
0] While the low-energy spectrum near the (magnetic) 
Bragg vectors looks conventional, the high-energy part 
of the excitation spectrum shows several anomalous fea- 
tures which both qualitatively and quantitatively differ 
strongly from LSWT: (i) At high energies the spectrum 
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is renormalized downwards with respect to LSWT. (ii) 
The excitation spectrum is very flat in the region halfway 
between the origin and a Bragg vector, with the mid- 
point renormalized downwards with respect to LSWT by 
~ 40%. (hi) There are local, roton-like, minima at the 
midpoints of the Brillouin zone edges, whose energy is 
~ 44% lower than the LSWT prediction. These strong 
downward renormalizations for the TLM should be con- 
trasted with the SLM, for which quantum fluctuations 
always renormalize the LSWT spectrum upwards and by 
amounts never exceeding 20%. 

Two-spinon interpretation. — We will argue that these 
results for the TLM are suggestive of spinons in the 
model. Our basic hypothesis is that an "uncorrelated" 
RPA-like calculation for the TLM, analogous to that 
discussed in Ref. for the SLM, should produce a 
spectrum similar to the LSWT result, but it will be 
modified by correlations [2l| . In particular, repulsion 
from the two-spinon (particle-hole) continuum can lower 
the magnon energy, especially at wavevectors where this 
continuum has minima. These minima should occur at 
{Ki — Kj)/2 corresponding to the creation of minimum- 
energy particle-hole excitations with particle and hole 
wavevectors Ki/2 and Kj/2, respectively, the locations 
of which are shown in Fig. 01 (the Ki are magnetic Bragg 
vectors). This equals (see Fig. 0J) a 'roton' wavevector 
when Ki and Kj differ by a 2ir/3 rotation around the 
origin, and a wavevector at a spinon minimum if Ki and 
Kj differ by a ir/3 rotation. Fig. ^ shows that at both 
types of wavevectors the excitation spectrum is strongly 
renormalized downwards with respect to the LSWT. At 
the former (latter) type of wavevector, the LSWT disper- 
sion is flat (peaked), which upon renormalization leads 
to a dip (flat region) in the true spectrum. Thus we at- 
tribute these deviations to the existence of a two-spinon 
continuum. 

For the SLM, the spectral weight of the magnon peak 
at (n, 0) is considerably smaller than at (tt/2,tt/2) (60% 
vs 85%) 0, and the magnon energy deviates much more 
from SWT at (tt,0) than at (7r/2,7r/2). This suggests 
quite generally that the relative weight of the magnon 
peak decreases with increasing deviation between the 
true magnon energy and the LSWT prediction. There- 
fore one might expect the contribution of the two-spinon 
continuum to be considerably larger for the TLM than 
for the SLM. 

As for the spinons proposed for the SLM 0, 0, 0] , 
the locations of the minima in the spinon dispersion re- 
flect a cJ-wave character of the underlying RVB pairing 
correlations. A <i-wave RVB state of this type, but with- 
out long-range order, was discussed for the TLM in Ref. 
I22I Its energy was however notably higher than that 
of the ordered ground state, and attempts to modify 
this RVB state to incorporate long-range order were not 
successful. A mean-field RVB state for the TLM with 
bosonic spinons, whose dispersion has minima close to 



Ki/2, was considered in Refs. 13 again without 
long-range order. In light of our spinon hypothesis for 
the TLM it would clearly be of interest to revisit these 
problems. 
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FIG. 4: Brillouin zone diagram for the TLM. Bragg vectors 
(filled squares), local roton minima (filled circles) and mid- 
points of flat dispersion regions (stars) for the spin-1 disper- 
sion are shown, as well as the global minima of the proposed 
5=1/2 spinon dispersion (open squares). Note that the lat- 
ter two sets of wavevectors coincide. The excitation spectra 
in Figs. and are plotted along the path ABCOAD. 

Explanation of finite-temperature anomalies. — The ex- 
istence of 'roton' minima and their description in terms 
of pairs of spinons provide a possible explanation for 
the sharply different temperature dependent properties 
of the SLM and TLM. For the SLM the temperature 
dependence of the correlation length is consistent with 
a NLSM description in the renormalized classical (RC) 
regime over a considerable temperature range |25| . Even 
though the ground state moment and spin-stiffness are 
comparable in the two models, for the TLM the correla- 
tion length was found to be orders of magnitude smaller 
at T = J/4 with a spin stiffness decreasing with decreas- 
ing temperature (and longer length scales) Q, inconsis- 
tent with the NLSM description in the RC regime [26| . 
We suggest that these differences are due to the fact that 
(see Figs. ^andOU the spinon gap E s (which is half the 
roton energy) is four times smaller for the TLM than for 
the SLM (0.28J versus 1.1J). Substantial thermal ex- 
citation of spinons for temperatures comparable to E s 
will make a significant contribution to the entropy and 
reduce the spin stiffness. Following an argument by Ng 
|27| . we expect that thermal excitation of spinons will 
cause a NLSM description to break down when T ~ E s . 
The results of the high-temperature expansions for both 
models are consistent with this estimate 

J\/J2 = 3 and Cs^CuCl^. — The ratio J1/J2 = 3, closer 
to the decoupled chains limit, is relevant for CS2CUCI4, 
which has an extremely rich excitation spectrum [jj] with 
well-defined spin-waves at low energies below the Neel 
temperature T^v, and a continuum, strongly reminiscent 
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of the two-spinon continuum in one-dimensional (ID) an- 
tiferromagnets, which persists well above T/v- Compared 
to LSWT the series dispersion (Fig. |3J is enhanced by 
~ 53% along the chains (close to the value for ID chains) 
and reduced by ~ 50% perpendicular to the chains. (In 
contrast, higher order SWT gives only weak enhance- 
ment (~ 13%) over LSWT along the chains.) Thus quan- 
tum fluctuations make the system appear much more ID. 
Overall, the series dispersion agrees well with the exper- 
imental dispersion for CS2CUCI4, also shown in Fig. 03 
whose enhancement and reduction factors (derived using 
a slightly different J\j Ji ratio; see caption) are ~ 63% 
and ~ 17%. The large difference between the theoretical 
and experimental reduction factors perpendicular to the 
chains may be due to the Dzyaloshinski-Moriya interac- 
tion in CS2CUCI4 (not included in the series calculations) 
which may make the system less ID (it has the same path 
as J2 and a coupling constant ~ O.I5J2) |28| . 
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FIG. 5: Excitation spectrum for J1/J2 = 3 (solid points 
from series) along the path ABCOAD in Fig. 0] compared to 
experimental dispersion in CS2CUCI4 (squares from Ref. Q, 
dashed line is experimental fit) where the exchange ratio is 
similar, J2/J1 = 0.34(3) and J 2 = 0.128(5) meV. Compared 
to LSWT for J1/J2 — S (full line) these spectra are enhanced 
along the Ji bonds and decreased perpendicular to them. 
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